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Abstract 

A new final state of gravitational collapse is proposed. By extending the concept of Bose- 
Einstein condensation to gravitational systems, a cold, dark, compact object with an interior de 
Sitter condensate p v = —p v and an exterior Schwarzschild geometry of arbitrary total mass M is 
constructed. These are separated by a shell with a small but finite proper thickness t of fluid with 
eq. of state p = +p, replacing both the Schwarzschild and de Sitter classical horizons. The new 
solution has no singularities, no event horizons, and a global time. Its entropy is maximized under 
small fluctuations and is given by the standard hydrodynamic entropy of the thin shell, which is of 
order k B £Mc/h, instead of the Bekenstein-Hawking entropy formula, S BH = Airk B GM 2 /he. Hence 
unlike black holes, the new solution is thermodynamically stable and has no information paradox. 
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Introduction. Cold superdense stars with a mass above some critical value undergo 
rapid gravitational collapse. Due to the impossibility of halting this collapse by any known 
equation of state for high density matter, a kind of consensus has developed that a collapsing 
star must inevitably arrive in a finite proper time at a singular condition, called a black hole. 

The characteristic feature of a black hole is its event horizon, the null surface of finite 
area at which outwardly directed light rays hover indefinitely. For simplicity, consider an 
uncharged, non-rotating Schwarzschild black hole with the static, spherically symmetric line 
element, 

dv^ 

ds 2 = -f(r) dt 2 + — + r 2 (d9 2 + sin 2 9 d^ 2 ) . (1) 
n{r) ' 

The functions fir) and h(r) are equal in this case, and 

f s (r) = h s (r) = l- r f, r s = 2 -^. (2) 

At the event horizon, r = r s , the metric (0) becomes singular. Since local curvature invari- 
ants remain regular at r = r s , a test particle falling through the horizon experiences nothing 
catastrophic there (if M is large enough), and it is possible to find regular coordinates which 
analytically extend the exterior Schwarzschild geometry through the event horizon into the 
interior region. 

It is important to recognize that this mathematical procedure of analytic continuation 
through a null hypersurface involves a physical assumption, namely that the stress-energy 
tensor is vanishing there. Even in the classical theory, the hyperbolic character of the 
Einstein equations allows generically for sources and discontinuities on the horizon which 
would violate this assumption. Whether such analytic continuation is mandatory, or even 
permissable in a more complete theory taking quantum effects into account, is still less 
certain. 

Non-analytic behavior is typical of quantum many-body systems at a phase transition. 
Quantum systems also exhibit macroscopic coherence effects which do not depend on local 
forces becoming large. Thus the fact that the tidal forces on classical test bodies falling 
through the event horizon are arbitrarily weak (proportional to M/r^ ~ 1/M 2 ) for an ar- 
bitrarily large black hole does not imply that quantum effects must be unimportant there. 
Electron waves restricted to the region outside an Aharonov-Bohm solenoid, where the elec- 
tromagnetic field strength vanishes and there are no classical forces whatsoever, nevertheless 
experience a shift in their interference fringes. Qualitatively new effects such as these arise 
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because quantum matter has extended wavelike properties, with many-body statistical corre- 
lations that can in no way be captured by consideration of pointlike test particles responding 
only to local forces. 

A photon with asymptotic frequency uj and energy hu far from the black hole has a 
local energy ftuj /~ 3 , which diverges at the event horizon. Unlike classical test particles, 
when h 7^ such extremely blue shifted photons are necessarily present in the vacuum as 
virtual quanta. Their effects upon the geometry depend upon the quantum state of the 
vacuum, defined by boundary conditions on the wave equation in a non-local way over all of 
space, and (XL 6 ) may be large at r = r„, notwithstanding the smallness of the local classical 

nrin 

curvature there |2J, |3J . Since the limits r — > r s and % — > do not commute, non-analytic 
behavior near the event horizon, quite different from that in the strictly classical (ft = 0) 
situation is possible in the quantum theory. 

The non-analytic nature of the classical limit H — > may be seen also in the thermody- 
namic analogy for the laws of black hole mechanics. Pursuing the analogy of these classical 
laws with thermodynamics, Bekenstein suggested that black hole event horizons carry an 
intrinsic entropy proportional to their area, Anr 2 jj, |f|. In order to have the units of 
entropy, the horizon area must be multiplied by a constant with units of c 3 k B /hG. When 
Hawking found that a flux of radiation could be emitted from a black hole with a well-defined 
temper-™, T. = t ,, W M fl g the thermodynamic analogy of area with entropy 
received some support, since the conservation of energy in the system can be written in the 
form, 

dE = T H dS BH = d {—j— ) , 0) 

suggestive of the first law of thermodynamics. 

The curious feature of relation (JHJ) is that H cancels out, and plays no dynamical role. 
The identification of S„„ with the entropy of the black hole is founded on the dynamics 

on n 

of purely classical relativity (i.e. % = 0), through the area law of refs. |9|, Hfl, . If this 
identification of rescaled classical area with entropy is to be valid in the quantum theory, 
then the limit h — > (with M fixed) which yields an arbitrarily low temperature, would 
assign to the black hole an arbitrarily large entropy, completely unlike the zero temperature 
limit of any other cold quantum system. 



Closely related to this paradoxical result is the fact, also pointed out by Hawking 12], 
that a temperature inversely proportional to M = E/c 2 implies a negative heat capacity: 



jS; = —(he 5 /8%GkB)T~ 2 < 0. However, a negative heat capacity is impossible for any 
system in stable equilibrium, since the heat capacity is proportional to the square of the 
energy fluctuations of the system. If this quantity is negative, the system cannot be in 
stable equilibrium at all, and the applicability of equilibrium thermodynamic relations is 
questionable. Attempts to evaluate the entropy of an uncharged black hole directly from 
statistical considerations (S = k B In Q) produce divergent results due to the unbounded 
number of wave modes infinitely close to the event horizon 14 1. The entropy of fields in 
a fixed Schwarzschild background would also be expected to scale linearly with the number 
N of independent fields, whereas S BH is independent of N. 

Ignoring these myriad difficulties and nevertheless interpreting (J3J) literally as a thermo- 
dynamic relation implies that a black hole has an enormous entropy, S BH ~ 10 77 k B (M/M Q ) 2 , 
far in excess of a typical stellar progenitor with a comparable mass. The associated 'infor- 
mation paradox' and implied violation of unitarity has been characterized as so serious as 
to require an alteration in the principles of quantum mechanics [15I . 

This paradoxical state of affairs arising from an area law originally derived in a strictly 
classical framework, together with the cancellation of h in . suggest that Q may not be 
a generally valid quantum relation at all, but only the (improper) classical limit of such a 
relation. The cancellation of h here is reminiscient of a similar cancellation in the energy 
density of modes of the radiation field in thermodynamic equilibrium, hu n(u) u 2 duj — > 
k B Tuj 2 duj in the Rayleigh- Jeans limit of very low frequencies, %uo <C k B T '. Improper extension 
of this low energy relation from Maxwell theory into the quantum high frequency regime leads 
to an ultraviolet catastrophe, similar to that encountered in the counting of wave modes near 
an event horizon. Conversely, treating the atoms in a solid as classical point particles leads 
to equipartition of energy and an incorrect prediction (the Dulong-Petit law) of constant 
specific heat for crystals at low temperatures. Both of these difficulties precipitated and 
were resolved by the rise of a new quantum theory of matter and radiation. 

In the illustrative case of the Einstein single frequency crystal, the difference between 
the heat capacity C V (T) at the temperature T and its Dulong-Petit high temperature limit, 
C v (oo) is AC V (T) = C V (T) — C v (oo), which has the leading high temperature behavior 
AC V (T) oc — T~ 2 , exactly the same temperature dependence as the negative Bekenstein- 
Hawking heat capacity of the Schwarzschild black hole. When the term C v (oo) arising from 
the discrete atomistic degrees of freedom is restored, the true heat capacity C V (T) of the 



4 



crystal is positive. Application of the Einstein energy fluctuation formula to a black hole 
shows that the Bekenstein-Hawking heat capacity is identical to what one would obtain for a 
large number of weakly interacting massive bosons, each with a mass of order of the Planck 
mass, close to the Dulong-Petit limit, but where the term analogous to C v (oo) has been 
dropped 

A related consideration suggesting that the Hawking emission depends on an improper 
ultraviolet extrapolation of classical waves into the quantum regime comes from examining 
the origin of these waves in the past. Although the Hawking temperature T H is very small for 
a large black hole, simple kinematic considerations show that the Hawking modes observed at 
a time t long after the formation of the hole originate as incoming vacuum modes at a greatly 
blue shifted frequency, uj in ~ (k B T H /h) exp(ct/2r 5 ). In other words, the calculation of the 
Hawking flux at late times assumes that the local properties of the fixed spacetime geometry 
are known with arbitrarily high precision, even at exponentially sub-Planck length and time 
scales, where one would normally question the semi-classical approximation, and indeed the 
existence of any well-defined metric at all. The contribution to the local stress-energy tensor, 
{T a b ) of these highly blue shifted modes diverges at the horizon, which is another form of 
the ultraviolet catastrophe. Only an exact balance with time-reversed highly blue shifted 
ingoing modes can cancel this divergence in the precisely tuned Hartle-Hawking thermal 
state [l8|. This state is in any case unstable due to its negative heat capacity. In any 
other thermal state, with T ^ T H , such as would be expected to arise from fluctuations, 
the divergence at the horizon is not cancelled, and one must expect a substantial quantum 
backreaction on the geometry there, which invalidates the basic assumption of an arbitrarily 
accurately known fixed classical Schwarzschild background, arbitrarily close tor = r s . 

Note that since it transforms as a tensor under coordinate transformations, a large lo- 
cal (T a b ) near r = r s is perfectly consistent with the Equivalence Principle, except in its 
strongest possible form which would admit no non-local effects of any kind, including those 
of macroscopic coherence and entanglement, known to exist in both relativistic and non- 
relativistic quantum many-body systems. 

In earlier models of backreaction the black hole was immersed in a Hawking radiation 
atmosphere, with an effective equation of state, p = up. It was found that due to the blue 
shift effect, the backreaction of such an atmosphere on the metric near r = r s is enormous, 
with the interior region quite different from the vacuum Schwarzschild solution, and with a 
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large entropy of order of S BH from the fluid alone [l^. In fact, S = Aj^S BH1 becoming 
equal to the Bekenstein-Hawking entropy for k — 1. Aside for accounting for the entropy 
S BH from purely standard hydrodynamical relations, this result suggests that the maximally 
stiff equation of state consistent with the causal limit may play a role in the quantum theory 
of fully collapsed objects. Despite such very suggestive features indicating the importance 
of backreaction on the geometry, these models cannot be viewed as a satisfactory solution to 
the final state of the collapse problem, since they involve huge (Planckian) energy densities 
near r s , as well as a negative mass singularity at r = 0. The negative mass singularity 
arises because a repulsive core is necessary to counteract the self-attractive gravitation of 
the dense relativistic fluid with positive energy. 



Recent 
made 



y 



another proposal for incorporating quantum effects near the horizon has been 
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231 ] , with a critical surface of a quantum phase transition replacing the classical 



horizon, and the interior replaced by a region with eq. of state, p = — p < 0. Equivalent to 
a positive cosmological term in Einstein's eqs^this eq. of state was first proposed for the 
final state of gra vitational collapse by Gliner [24] and considered in a cosmological context 
by Sakharov |2£|. It violates the strong energy condition p + 3p > used in proving the 
classical singularity theorems. As is now well known because of the observations of distant 
supernovae implying an accelerating universe, the geodesic worldlines of test particles in a 
region of spacetime where p + 3p < diverge from each other, mimicking the effects of a 
repulsive gravitational potential. Thus such a region free of any singularities in the interior 
region r < r s can replace the unphysical negative mass singularity encountered in the pure 
K — 1 fluid model. 

Based upon these considerations, in this paper we show that a consistent static solution 
of Einstein's equations can be constructed, with the critical surface of refs. 21, 2^| replaced 
by a thin shell of ultra-relativistic fluid with eq. of state p = p. Because of its vacuum energy 
interior with p = —p the new solution is stable and free of all singularities. Its entropy is a 
local maximum of the hydrodynamic entropy, whose modest value given by (j27)l - (j28|l below 
is easily attainable in a physical collapse process from a stellar progenitor with a comparable 
mass. 

The model we arrive at is that of a low temperature condensate of weakly interacting 
massive bosons trapped in a self-consistently generated cavity, whose boundary layer can 
be described by a thin shell. The assumption required for a solution of this kind to exist is 
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that gravity, i.e. spacetime itself, must undergo a quantum vacuum rearrangement phase 
transition in the vicinity of r = r s . In this region quantum zero-point fluctuations dominate 
the stress-energy tensor, and become large enough to influence the geometry, regardless of the 
composition of the matter undergoing the gravitational collapse. As the causal limit p = p 
is reached, the interior spacetime becomes unstable to the formation of a gravitational Bose- 
Einstein condensate (GBEC), described by a non-zero macroscopic order parameter in the 
effective low energy description. Since a condensate is a single macroscopic quantum state 
with zero entropy, a model of a cold condensate repulsive core as the stable, non-singular 
endpoint of gravitational collapse provides a resolution of the information paradox, which 
is completely consistent with quantum principles j^J. The interior and exterior regions 
are separated by a thin surface layer near r = r s where the vacuum condensate disorders. 
Any entropy in the configuration can reside only in the excitations of this boundary layer. 
A suggestion for the effective theory incorporating the effects of quantum anomalies that 
describes this fully disordered phase, where the role of the order parameter is played by 
the conformal part of the metric has been presented elsewhere j^J. For a recent review of 
investigations of other non-singular quasi-black-hole (QBH) models see 

The Vacuum Condensate Model. In an effective mean field treatment for a perfect 
fluid at rest in the coordinates (0), any static, spherically symmetric collapsed object must 
satisfy the Einstein eqs. (in units where c = 1), 

-G\ = ~ [r (1 - h)] = -SitGT\ = SttGp, (4) 

G r r = — £ d f + ^(h-l) = ^GT r r = MGp, (5) 
rj ar r 2 

together with the conservation eq., 

V„T" r = ^ + ^f + 2^ = 0, (6) 
ar 2/ ar r 

which ensures that the other components of the Einstein eqs. are satisfied. In the general 
spherically symmetric situation the tangential pressure, pj_ = T e e = T"^ is not necessarily 
equal to the radial normal pressure p = T r r . However, for purposes of developing the simplest 
possibility first, we restrict ourselves in this paper to the isotropic case where p± = p. In 
that case, we have three first order eqs. for four unknown functions of r, viz. f, h, p, and p. 
The system becomes closed when an eq. of state for the fluid, relating p and p is specified. 



Because of the considerations of the Introduction we allow for three different regions with 
the three different eqs. of state, 

I. Interior : < r < r\ , p = —p , 

II. Thin Shell : r\ < r < r 2 , p — +p , (7) 

III. Exterior : r 2 < r , p = p = . 

In the interior region p = —p is a constant from (jUJ). Let us call this constant p v = 
3Hq/8ttG. If we require that the origin is free of any mass singularity then the interior is 
determined to be a region of de Sitter spacetime in static coordinates, i.e. 

I. f{r) = Ch(r) = C(l -H 2 r 2 ), < r < r x , (8) 

where C is an arbitrary constant, corresponding to the freedom to redefine the interior time 
coordinate. 

The unique solution in the exterior vacuum region which approaches flat spacetime as 
r — > oo is a region of Schwarzschild spacetime (J2J), viz. 

HI. f( r ) = h{r) = l--^—, r 2 <r. (9) 

The integration constant M is the total mass of the object. 

The only non-vacuum region is region II. Let us define the dimensionless variable w by 
w = 8nGr 2 p, so that eqs. (JU)-© with p = p may be recast in the form, 

dr dh 



r 1 — w — h 
dh 1 — w — h dw 



h 1 + w — 3h w 



(11) 



together with pf oc wf /r 2 a constant. Eq. (fTU|) is equivalent to the definition of the 
(rescaled) Tolman mass function by h = 1 — 2m(r)/r and dm(r) = AirG pr 2 dr = wdr/2 
within the shell. Eq. can be solved only numerically in general. However, it is possible 
to obtain an analytic solution in the thin shell limit, < h <C 1, for in this limit we can set 
h to zero on the right side of to leading order, and integrate it immediately to obtain 

2m (1 + w^ 2 



— ~e- -<1, (12) 

r w 



in region II, where e is an integration constant. Because of the condition h <C 1, we require 
e<l, with w of order unity. Making use of eqs. ([TU J) -(fTT |) and (fT2j) we have 

±*-edw£±^-. (13) 

Because of the approximation e«l, the radius r hardly changes within region II, and dr is of 
order edw. The final unknown function / is given by / = (r/r 1 ) 2 (w 1 /if)/(r 1 ) ~ (wi/w)f(ri) 
for small e, showing that / is also of order e everywhere within region II and its boundaries. 

At each of the two interfaces at r = r\ and r = r 2 the induced three dimensional metric 
must be continuous. Hence r and f(r) are continuous at the interfaces, and 

ft \ Wl ft \ CWl (^ C/2 2x , 2GM /-.^ 

/(r 2 ) ~ — /(ri) = l-# n =1 . 14 

To leading order in e -C 1 this relation implies that 

n ~ ^- ~ 2GM ~ r 2 . (15) 

Thus the interfaces describing the phase boundaries at ri and r 2 are very close to the classical 
event horizons of the de Sitter interior and the Schwarzschild exterior. 

The significance of < e -C 1 is that both / and h are of order e in region II, but 
are nowhere vanishing. Hence there is no event horizon, and t is a global time. A photon 
experiences a very large, 0(e~^) but finite blue shift in falling into the shell from infinity. 
The proper thickness of the shell between these interface boundaries is 

T2 1 1 Z^ 1 3 1 

c?r/i _ 2~r s e2 / dw w~z ~ e^r s , (16) 
n J w 2 

and very small for e — > 0. Because of (fTB^) the constant vacuum energy density in the interior 
is just the total mass M divided by the volume, i.e. p v ~ 3M/47rr^, to leading order in e. 
The energy within the shell itself, 

f 2 f Wl dw 

En = 4vr / pr 2 dr ~ eM I — (1 + w) ~ eM (17) 

is extremely small. 

We can estimate the size of e and i by consideration of the expectation value of the 
quantum stress tensor in the static exterior Schwarzschild spacetime. In the static vacuum 
state corresponding to no incoming or outgoing quanta at large distances from the object, 
i.e. the Boulware vacuum , the stress tensor near r = r s is the negative of the stress 



tensor of massless radiation at the blue shifted temperature, T\ oc —T H j \J f{r) and diverges 
as Ti oc ~ /~ 2 (?") as r — * r s . The location of the outer interface occurs at an r where this 
local stress-energy oc M~ 4 e -2 , becomes large enough to affect the classical Schwarzschild 



curvature ~ M 2 , i.e. when 



P l 



M 



10 



-38 



Mr. 







M 



(18) 



where Mpi is the Planck mass \Jhc/G ~ 2 x 10~ 5 gm. Thus e is indeed very small for 
a stellar mass object, justifying the approximation a posteriori. With this semi-classical 
estimate for e we find 



PI 1 s — 



3 x 10 



-14 



M 



cm. 



(19) 



Although still microscopic, the thickness of the shell is very much larger than the Planck 
scale Lpi ~ 2 x 10~ 33 cm. The energy density and pressure in the shell are of order M~ 2 
and far below Planckian for M 3> Mpi, so that the geometry can be described reliably by 
Einstein's equations in both regions I and II. 

Although f(r) is continuous across the interfaces at r\ and r2, the discontinuity in the 
eqs. of state does lead to discontinuities in h(r) and the first derivative of f(r) in general. 
Defining the outwardly directed unit normal vector to the interfaces, n b = 5 r b ^Jh{r), and 
the extrinsic curvature K a b = V a n b , the Israel junction conditions determine the surface 
stress energy rj and surface tension a on the interfaces to be given by the discontinuities in 
the extrinsic curvature through 



~2f~dr 



4irG(r) - 2a) 



[K e e ] = [Ka 



-AttGt] . 



(20) 



(21) 



Since h and its discontinuities are of order e, the energy density in the surfaces rj ~ ea, while 
the surface tensions are of order e~2. The simplest possibility for matching the regions is to 
require that the surface energy densities on each interface vanish. From ()21)1 this condition 
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implies that h(r) is also continuous across the interfaces, which yields the relations, 



From (|13|) dw/dr < 0, so that w 2 < w\ and C < 1. In this case of vanishing surface energies 
rj = the surface tensions are determined by (J20|) - (}2T|) to be 



to leading order in e at r± and r 2 respectively. The negative surface tension at the inner 
interface is equivalent to a positive tangential pressure, which implies an outwardly directed 
force on the thin shell from the repulsive vacuum within. The positive surface tension on 
the outer interfacial boundary coresponds to the more familiar case of an inwardly directed 
force exerted on the thin shell from without. 

The entropy of the configuration may be obtained from the Gibbs relation, p+p = sT+np, 
if the chemical potential p is known in each region. In the interior region I, p + p — 
and the excitations are the usual transverse gravitational waves of the Einstein theory in 
de Sitter space. Hence the chemical potential p may be taken to vanish and the interior 
has zero entropy density s = 0, consistent with a single macroscopic condensate state, 
S = k B lnQ = for Q = 1. In region II there are several alternatives depending upon 
the nature of the fundamental excitations there. The p = p eq. of state may come from 
thermal excitations with negligible p or it may come from a conserved number density n of 
gravitational quanta at zero temperature. Let us consider the limiting case of vanishing p 
first. 

If the chemical potential can be neglected in region II, then the entropy of the shell is 
obtained from the eq. of state, p = p = (a 2 /8irG)(k B T/h) 2 . The T 2 temperature dependence 
follows from the Gibbs relation with p = 0, together with the local form of the first law 
dp = Tds. The Newtonian constant G has been introduced for dimensional reasons and a is 
a dimensionless constant. Using the Gibbs relation again the local specific entropy density 




(23) 



(24) 



(22) 




(26) 



(25) 
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s(r) = a 2 k 2 B T{r) /Anh 2 G = a{k B /Ti){p/2TiG)^ for local temperature T(r). Converting to our 
previous variable w, we find s = (ak B / AnhGr) and the entropy of the fluid within the 
shell is 

S = An P sr 2 drh-^ a ^ie^ In V (27) 
to leading order in e. Using (JTfiJl and (JUJ), this is 
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S-^K^-^ 57 ak B (^-j «S flH . (28) 

The maximum entropy of the shell and therefore of the entire configuration is some 20 orders 
of magnitude smaller than the Bekenstein-Hawking entropy for a solar mass object, and of 
the same order of magnitude as a typical progenitor of a few solar masses. The scaling of 
(I28|) with M§ is also the same as that for supermassive stars with M > 100 M , whose 
pressure is dominated by radiation pressure 30]. Thus the formation of the GBEC star 
from either a solar mass or supermassive stellar progenitor does not require an enormous 
generation or removal of entropy, and there is no information paradox. 

Because of the absence of an event horizon, the GBEC star does not emit Hawking 
radiation. Since w is of order unity in the shell while r ~ r s , the local temperature of 
the fluid within the shell is of order T H ~ h/k B GM. The strongly redshifted temperature 
observed at infinity is of order yfeT H , which is very small indeed. Hence the rate of any 
thermal emission from the shell is negligible. 

If we do allow for a positive chemical potential within the shell, \x > 0, then the temper- 
ature and entropy estimates just given become upper bounds, and it is possible to approach 
a zero temperature ground state with zero entropy. This non-singular final state of gravita- 
tional collapse is a cold, completely dark object sustained against any further collapse solely 
by quantum zero-point pressure. 

Stability. In order to be a physically realizable endpoint of gravitational collapse, any 
quasi-black hole candidate must be stable Q| • Since only the region II is non- vacuum, with a 
'normal' fluid and a positive heat capacity, it is clear that the solution is thermodynamically 
stable. The most direct way to demonstrate this stability is to work in the microcanonical 
ensemble (in the case of zero chemical potential) with fixed total M, and show that the 
entropy functional, 

n ak R n , ( dm\^ f 2m(r)\^ . . 

S=t£ rdr 2— 1-—^ , (29) 



hG J ri \ dr J 

12 



is maximized under all variations of m(r) in region II with the endpoints (ri,r2), or equiv- 
alently {w^Wq) fixed. 

The first variation of this functional with the endpoints r\ and r 2 fixed vanishes, i.e. 5S = 
by the Einstein eqs. Q-© for a static, spheric ally symmetric star. Thus any solution of 
eqs. P|)-(|fi|) is guaranteed to be an extremum of S 31] . This is also consistent with regarding 
Einstein's eqs. as a form of hydrodynamics, strictly valid only for the long wavelength, 
gap less excitations in gravity. In the context of a hydrodynamic treatment, thermodynamic 
stability is also a necessary and sufficient condition for the dynamical stability of a static, 
spherically symmetric solution of Einstein's equations j^ . 

The second variation of ()29jl is 



{- 


d(5m) 




dr 



when evaluated on the solution. Associated with this quadratic form in 5m is a second order 
linear differential operator £ of the Sturm-Liouville type, viz. 
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^ _ d J / c?m\ 2 I _|_ 2 / dm\ 2 ( ^ ^ <im\ 
c?r | \ dr ) dr J r \ dr J \ dr J 

This operator possesses two solutions satisfying Cx = 0, obtained by variation of the 
classical solution, m(r; r 1; r 2 ) with respect to the parameters (r 1? r 2 ). Since these correspond 
to varying the positions of the interfaces, x does not vanish at (r 1 ,r 2 ) and neither function 
is a true zero mode. For example, it is easily verified that one solution is Xo = 1 — w -> from 
which the second linearly independent solution (1 — w) In w + 4 may be obtained. For any 
linear combination of these we may set 5m = x i>, where ip does vanish at the endpoints 
and insert this into the second variation ()30|) . Integrating by parts, using the vanishing of 
5m at the endpoints and Cx = gives 



ri 



Thus the entropy of the solution is maximized with respect to radial variations that vanish 
at the endpoints, i.e. those that do not vary the positions of the interfaces. Perturbations 
of the fluid in region II which are not radially symmetric decrease the entropy even further 
than (}3*2*|) . which demonstrates that the solution is stable to all small perturbations keeping 
the endpoints fixed. 
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Allowing for endpoint variations as well requires the inclusion of the vacuum stress, (T a b ) 
in the vicinity of the interfaces, which fixed e by the estimate (j!8|) . It is clear that the 
vacuum (T a b ) must be included in a more complete model for another reason. The general 
stress-energy in a spherically symmetric, static spacetime has three components, namely p,p 
and p±. We have set p± = p and restricted ourselves to only two isotropic eqs. of state 
p = —p and p = +p only for simplicity, to illustrate the general features of a non-singular 
solution to the gravitational collapse problem in a concrete example. For any static solution, 
we must expect also the stress-energy tensor of vacuum polarization in the Boulware vacuum 
to contribute. This stress-energy satisfies p = p/3 < near the horizon The addition of 
such a negative pressure eq. of state in the thin outer edge obviates the need for the positive 
pressure discontinuity from negative pressure inside to positive pressure outside. Hence a 
completely smooth matching of h and df/dr is possible and the surface tensions (|25 p -([26 |) 
can be made to vanish identically. A full analysis of dynamical stability without restriction 
on the interface boundaries will be possible in the framework of a more detailed model which 
leads to these vacuum stresses in the boundary layer. Such an investigation can be carried 
out without reference to thermodynamics or entropy and would apply then even in the case 
of a configuration at absolute zero. 

Conclusions. A compact, non-singular solution of Einstein's eqs. has been presented 
here as a possible stable alternative to black holes for the endpoint of gravitational collapse. 
Realizing this alternative requires that a quantum gravitational vacuum phase transition 
intervene before the classical event horizon can form. Since the entropy of these objects is 
of the same order as that of a typical stellar progenitor, even for M > 1OOM , there is no 
entropy paradox and no significant entropy shedding needed to produce a cold gravitational 
vacuum or 'grava(c)star' remnant. 

Since the exterior spacetime is Schwarzschild until distances of order of the diameter of 
an atomic nucleus from r = r s , a gravastar cannot be distinguished from a black hole by 
present observations of X-ray bursts 32]. However, the shell with its maximally stiff eq. 
of state p = p, where the speed of sound is equal to the speed of light, could be expected 
to produce explosive outgoing shock fronts in the process of formation. Active dynamics 
of the shell may produce other effects that would distinguish gravastars from black holes 
observationally, possibly providing a more efficient particle accelerator and central engine 
for energetic astrophysical sources. The spectrum of gravitational radiation from a gravastar 
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should bear the imprint of its fundamental frequencies of vibration, and hence also be quite 
different from a classical black hole. 

Quantitative predictions of such astrophysical signatures will require an investigation of 
several of the assumptions, and extension of the simple model presented in this paper in 
several directions. Although the eq. of state p = p is strongly suggested both by the limit of 
causality characteristic of a relativistic phase transition, and by the correspondence of the 
fluid entropy with S BH when the inner GBEC region is shrunk to zero, this eq. of state has 
been assumed here, not derived from first principles. Knowledge of the effective excitations 
in the shell is necessary to determine the chemical potential /i, and whether the entropy 
estimate (|28|) is accurate, or more properly to be regarded as an upper bound on the entropy 
of a GBEC star. The neglect of the p = p/3 < vacuum polarization in our model leads to 
some freedom in matching at the two interfaces and the surface tensions (l25j) - (f2T)j) . which may 
be different in detail or not present at all in a more complete treatment. A full analysis of 
the dynamical stability of the object, including the motion of the interfaces or the boundary 
layer(s) which replace them requires at least a consistent mean field description of quantum 
effects in this transition region. Although general theoretical considerations indicate that 
non-local quantum effects may be present in the vicinity of classical event horizons, a detailed 
discussion of how these effects can alter the classical picture of gravitational collapse to a 
black hole has not been attempted in this paper. Lastly, distinguishing the signatures of 
gravastars from classical black holes in realistic astrophysical environments, such as in the 
presence of nearby masses or accretion disks will depend on the details of the dynamical 
surface modes, as well as the extension of the spherically symmetric static model presented 
here to include rotation and magnetic fields. 

One may regard the model presented in this paper as a proof of principle, the simplest 
example of a physical alternative to the formation of a classical black hole, consistent with 
quantum principles, which is free of any interior singularity or information paradox. Ad- 
ditional theoretical and observational effort will be required to establish the cold, dark, 
compact objects proposed in this paper as the stable final states of gravitational collapse. 

Finally let us note that the interior de Sitter region with p = —p may be interpreted 
also as a cosmological spacetime, with the horizon of the expanding universe replaced by 
a quantum phase interface. The possibility that the value of the vacuum energy density 
in the effective low energy theory can depend dynamically on the state of a gravitational 



15 



condensate may provide a new paradigm for cosmological dark energy in the universe. The 
proposal that other parameters in the standard model of particle physics may depend on 
the vacuum energy density within a gravastar has been discussed by Bjorken 33]. 
Research of P. O. M. supported in part by NSF grant 0140377. 
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